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Abstract: The correlation function oi aV — A current with aV-\-A current is discussed 
within the framework of QCD in the Hmit of a large number of colours Nc- Applications 
to the evaluation of chiral condensates of dimension six and higher, as well as to the 
matrix elements of the Qj and Qs electroweak penguin operators are discussed. A critical 
comparison with previous determinations of the same parameters has also been made. 
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1. Introduction 

We shall be concerned with the correlation function of a left-handed current with a right- 
handed current 

Lt'ix) =u{x)j^^^^^d{x) and R'^ (Q) = d{0)Y ^-^u{0) , (1.1) 

in QCD and in the chiral limit where the light quarks are massless. In this limit, the 
correlation function in question depends only on one invariant amplitude liiniQ^) the 
euclidean momentum squared = —q^, with q the momentum flowing through the two- 
point function (see Fig. 1): 




Fig. 1 Feynman diagram representing the Ulr correlation function in Eq. in the 

large-Nc limit. The solid lines are light quarks propagating in a gluonic background of 
large-Nc planar diagrams. 

The interest on the function IIlr{Q'^) is twofold: on the one hand, as pointed out in 
ref. |l|, liiRiQ^) in the chiral limit is an order parameter of spontaneous chiral symmetry 
breaking (SxSB) at all values of Q'^; on the other hand, part of the absorptive part of this 
function i.e., ^lmIlLR{t) with t = > 0, is accessible to experiment via hadronic r-decays 
and e'^e~ annihilation into hadrons. Furthermore, it has also been shown ^, ^, ^, that 
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the same IIlb.{Q'^) function governs the underlying dynamics of the leading contributions 
to the weak matrix elements of the electroweak Penguin-like operators 



Qj = Q{sL7^dL) ^ egiqRj^qR) and Qs = -12 ^ eg{sLqR){qRdL) , (1-3) 

q=u,d,s q=u,d,s 

where eg denote quark electric charges in units of the electric charge and summation over 
quark color indices within brackets is understood. These considerations raise the liLRiQ"^) 
function to the rank of an excellent theoretical laboratory to test new ideas on the funda- 
mental subject of SxSB in QCD. 

Here, we shall be particularly concerned with the study of Hlr{Q^) iu the limit of a 
large number of colors Nc in QCD. Part of our motivation is to understand the discrepancies 
between various phenomenological analyses which have recently been made g 0, |ll|, |l| 



|l^, 0, |T5| using the experimental data on hadronic r-decays published by the ALEPH Q 
and OPAL Q collaborations at LEP. There are conflicting results for the chiral condensates 
which modulate the asymptotic behaviour of the Il-LRiQ^) function at large values, 
between those obtained in refs. |^, 10, 11, and those obtained in refs. |12, 13, 14]. We 



want to compare these results to those obtained in two successive approximations to the 
large-A'^c limit: the so-called minimal hadronic approximation (MHA) [|lO| consisting of 
a spectrum of a pion state, a vector state and an axial vector state; and the improved 
approximation where an extra higher vector state is added. 

In practice, we shall be working with the dimensionless complex function de- 
fined as 

Wlr[z] = -zHLRizM^p) , with Rez = ^, (1.4) 

and use the mass of the lowest massive state, the p(770 MeV), to normalize quantities 
with dimensions. In large- A'^c QCD the function VFlr[z] is a meromorphic function and, 
therefore, in full generality, it can be approximated by successive partial fractions of the 
type 

p ^ N 

Wlr[z] = AnY]^!— — tI1(2 + '7j), with pi = l and Pi ^ Pk for i ^ k ; (1.5) 

i=i y^ + Pi) j=i 

where P (and A^) get larger and larger, but finite. On the other hand, in QCD, the 
operator product expansion (OPE) of the two currents in Eq. (1^) fixes the large-Q^ fall 
off in l/Q^-powers of the invariant function HiRiQ^) to 

^lhn^n^^(Q2) = 5^^C2„+4(Q^/.^)(02„+4(^^))^^2)^ = 2 (^-6) 



Matching the leading asymptotic behaviour for large-z in Eq. ( |l.q ) to the one of the OPE 
in Eq. (|1.(]|), restricts the number of zeros A^ and the number of poles P in Eq. ( |1.E| ) to 
obey the constraint 

N-P = -2. (1.7) 
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In the chiral limit, all the vev's {02n+4) in Eq.(|1.6|) are order parameters of SxSB. In 
particular, the normalization factor A]\f in Eq. ( |1.5| ) is fixed by the residue of the leading 
term in the OPE, 

The case where = corresponds to the MHA already mentioned where, besides 
the Goldstone pole (the pion), which has been removed by the factor in the r.h.s. of 
Eq. (1.4), there are two poles (since P = 2 in this case): the lowest vector state and an 
axial-vector state [j^]. 

Strictly speaking, the (C'2ra+4) in the r.h.s. of Eq. (|1.6|) are dependent, because of 
the logQ^ dependence of the Wilson coefficients C2n+4(Q^, /U^) via the pQCD series in Us- 
This logQ^ dependence can only be reproduced if P ^ oo (and N cc) in Eq. (|1.5| ). In 



writing the large- A'^c approximation to the W^ijjz] function in Eq. (1.5), one is implicitly 
assuming an effective cancellation between the extra poles and zeros in the complex z- 
plane which lye beyond a disc of radius sq covering all the poles and zeros retained in 
that approximation. In the Minkowski axis, this is equivalent to assuming that, for all 
practible purposes, there is a cancellation between the ^-spectral function and the A- 
spectral function for t > sq, an assumption which is consistent with the fact that in pQCD 
and in the chiral limit these two spectral functions are identical. The same approximation, 
in the deep euclidean region, suggests that when comparing the asymptotic inverse powers 
of z in Eq. (|1.5| ) to those of the OPE in Eq. ( |1.6| ), the logQ^ dependence of the Wilson 
coefficients in C2n+4(Q^, /U^) be frozen at ~ sq. 



2. General Properties of IIlr{Q^) in Large-A^^ QCD 

The function I[lji{Q'^) is the simplest case of a class of Green's functions, which we call 
generalized two-point functions, and which obey rather remarkable short-distance ^ long- 
distance duality properties. Generalized two-point functions are two-point functions with a 
finite number of local operator insertions carrying zero momenta, (no insertions in the case 
of IIlii{Q'^)). These Green's functions, in the lavge-Nc limit, are meromorphic functions 
which can be approximated by polynomial ratios 

p ^ N 

W[z; pi,p2,---PP; 0-i,cr2---crAr] =^Ar JJ— rJI^^ + '^i)! (^•^) 

i=l + j=l 

and they are fully determined by their poles, their zeros and the overall normalization Ajsf. 
Using this representation, one can readily see that 



W 



111 111 1 

z' Pl P2 PP' CTi' a2 CTN 



P—N llPi Tirr 1 

z — W[z; pi,p2,--- Pp', (^i,o-2---<yN\, 

(2.2) 

relating the OPE expansion (powers of 1/z) to the chiral expansion (powers of z). In the 
case of the correlation function in Eq. ( [1.4| ) this is the generalization to an arbitrary 
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number of narrow states of the simple relation 



Wlr 



11 1 

Pi P2 



PiP2Wlr[z; pi,P2], for N = 0, (2.3) 



corresponding to the case of the MHA, where P = 2 and = 0. 

One can easily deduce the following set of properties, useful for phenomenological 
applications: 

1. The function at the origin 

WlrM = ^^ ^i^2.--cTjv ^ ^ ^ ^^^^ Wlr[0] = for = (MHA) . 

PlP2 ■■■PP PlP2 

(2.4) 

Recall that the value of Wl_r[0] is fixed by lowest order chiral perturbation theory; so, 
once An is known, the sign of Wlr[0] fixes the sign of the product of all the zeros ^. 
We therefore have that 



A^^^^^^^ = ^ ^ PF for Ar>l and = ^ for A = 0, (2.5) 

PIP2---PP M{, piP2 



with Fq the residue of the Goldstone pole (the pion) contributing to ^LpiQ"^)- 
In what follows, we shall often use the notation p2 = pa and p3 = py- 

2. The Wl/j[2;] function at infinity: condensates 

Starting from equation (|1.5D with the constraint N — P = —2 (which means that the 



two Weinberg sum rules are automatically satisfied) and decomposing it into partial 
fractions, we obtain 



p-2 



i=i ^ j=i k=i ^'^ 



with 



(2.7) 

where 5ik is the Kronecker symbol. This fixes the residues of the poles, which are 
physical couplings, in terms of the poles and zeros. There follows then that: 

p 

(02n+4) = 2(-l)"+i m2"+4 ""^pI ■ (2-8) 

fc=l 



^Notice that, if there are complex zeros, they have to appear in conjugate pairs of each other and 
therefore give a positive contribution to the product. 
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The Linear Constraint 

This is a very interesting constraint, which simply follows by expanding Eq. ( |1.5| ) to 
first non-trivial order in inverse powers of z 

It relates the difference between the sum of the positions of the zeros and the sum of 
the positions of the poles to the ratio of next -to-leading to leading vev's in the OPE. 
In other words, once we know the positions of the poles, the sum of the positions of 
the zeros, which is always real, is governed by the next-to-leading order term in the 
OPE. In the case corresponding to the MHA, where by definition there are no zeros, 
this constraint simply becomes: 



1 (O; 



implying that, in the MHA (Os) and (Oe) must have opposite signs. This is indeed 
what the phenomenological analyses in refs. 0, |l^ find, in contradistinction to 

the results in r efs. P, H, 11 which find the same sign for the first two condensates. 
In fact, in the case of the MHA, we can show that this alternance of sign for two 
consecutive condensates is a generic property. It follows from the relation 

(02n+4) =2(-irM2"+Mo^^ = 2(-irM2-+VF Vpi, for n>l. 

(2.11) 

Except for the (—1)" factor, all the quantities appearing in the r.h.s. are positive; 
which explains the alternance in sign with n. 

4. The Slope Constraint 

The value of the derivative of Wi:_R[2;] at z = is controlled by the 0{p'^) xPT low 
energy constant Liq [^] 

T^i^[0]=4Lio. (2.12) 

The relation to the poles and zeros of the I^i^ij[2;] -function in large- A'^c QCD is as 
follows 

( ^ 1 ^ l\ 
WU[0]=Pf\Y.--J2-\ ^ ^^1- (2.13) 



In the case A^ = 1 it reduces to the interesting relation 



''-Ml 



1/11 
-- 1 + — + — 

<y \ PA Pv 



(2.14) 



and in the case A = 0, corresponding to the MHA, to the well-known result |1£, 20| 



4Lio = f 1 + — ) ^ Lio = -^^ for pA = 2, (2.15) 
\ PaJ 8M^ 
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indicating that in the MHA, the slope of M/l_r[2;] at the origin has the opposite sign 
to the value of V^i,_r[0], in agreement with experiment. 



5. Dispersion Relations 

The function VFififz] and its corresponding spectral function ^lmIlLR{t) are related 
to each other by the dispersion relation 

O'^ f°° 1 1 

T^iR(Rez) = -^ / dt . -lmn^R(t). (2.16) 

Jo t + Q-^ — le-K 

This is the Hilbert transform of a spectral function, which in our approximation of 
large-A^c QCD has the general form 

-ImULRit) = -F^5{t) + V aiMfS{t - Mf) (2.17) 



1=1 



with the sum ordered in increasing values of the masses Mf. The residues ai are 
positive for the vector states and negative for the pion pole (the first term) and the 
axial states. 

Sometimes it is also convenient to consider the Laplace transform (Borel transform) 
of the spectral function 

1 r°° 1 

MLR{a) = Y72 dte-''' -ImULRit). (2.18) 

J^p Jo TT 

The two types of transforms are related by the fact that 



t + 

therefore. 



1 



(2.19) 



poo 

Wlr{z) = -zMl / dae-'''^>MLR{<j) . (2.20) 

^0 

In practice, one is often interested in observables which are moments of Wlr{Q^) 

(2.21) 

showing that we can view the observables 0^'^\ either as moments of Wlr{z), or as 
inverse moments of MLR{cr)- 

Some interesting properties of the Laplace Transform follow. 

• The Laplace transform M-LRic] obeys the differential equation 

|ap^ + ap_i^-p3j- + --- + ai^+ao|A^Li?H =0, (2.22) 
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with ap, ap-i, ... qq the coefficients of the polynomial 
p 



Y\_{z + Pi) = apz^ + ap-iz^ ^ H h oq ; 



(2.23) 



i=l 



I.e., 



p 

ap = l, ap_i = ^/>i, ••• ao = JJ/9i. (2.24) 

i=l i=l 

The polynomial Hilil-^ + Pi) is a stable polynomial., because all its roots are in 
the negative real axis (the Minkowski axis). 

The most general solution of the differential equation in ( |2.22| ) is of the form 



El 



+ aipie P'^" + a2P2e H \-apppe 



'Ppcr 



(2.25) 



with the Ui constants fixed by the successive boundary conditions obtained e.g., 
from the knowledge of the derivatives of the function Mlr[o'] at the origin. This 
leads to the system of equations first discussed in ref. ||l|: 



f p2 

^+ aipi + 02/02 H h appp = Mlr[Q] , 



aiPi + 02/02 "I 1" Oippp 



dMi 



(2.26) 



P+l , P+l , , P+l d^MLRlO] 

aiPi ^ + a2p2 +■■■ + appp^ = ■ 



The discriminant of this system of equations in the Oj is the Vandermonde 
determinant 



/II. 


■ 1 \ 


pi P2 ■ 


• PP 


P P 
\Pl P2 ■ 


• Pp) 



and 



d^Mr 



On the other hand the successive values of M lr [0] , ^"^^^ , • • • '^^^^ — (da)P~ 
are fixed by the OPE of the function Wlr[z]. In our case A^Lij[0] = 0; = 
and the first two equations in ( |2.26| ) are nothing but the well-known 1st and 
2nd Weinberg sum rules. 

• The positive moments of the Laplace transform, with the Goldstone singularity 
removed, correspond to coefficients of the chiral expansion 



daa^MiRia) = T{N + 1) / dt-^-lmULR{t) , for iV > , (2.27) 



1 1_ 



where 



-ImllLRit) = -ImULRit) + F^Sit) . 

TT TT 



(2.28) 
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6. Reconstruction of the Spectral Function 

It is useful to give the coefficients Qj of equation (|2.17] ) in terms of the poles pi and 



Inserting ( p^D into ( |2l6| ) we find 



Wlr[z\ = Pf - ^ OLkPk + ^ 



1 



k=l 



Z + Pk 



By comparing (|2.6| ) and ( ^.29 ), we then have 



- ^ OikPk = 0, 



k=l 



(2.29) 



(2.30) 



"fcPi = w^fc , with Y.k=i Wk = 0, 

where the two sum equations above are just the first and second Weinberg sum rules. 
Using a more conventional notation in terms of axial to vacuum and vector to vacuum 



couplings: ua = —f\ and ay = fy, we conclude from equations ( 2.30 ) and ( p.7[ ) that 



fvM^ 



-^0 TT Pi 

PV fj^ {Pi - PV + ^iv) 



P-2 

n 



(cTj - pv) 



(2.31) 



4n 



p-2 

n 



(cTj - pa) 



PA fJl {Pi - PA + SiA) fj^ 



(2.32) 



7. The Smoothness Assumption 

As already stated, the MHA corresponds to the case where = 0. In this ap- 
proximation, the function VF[z], or its equivalent M{a), is a monotonous function 
of the euclidean variable: < < cxd or, equivalently, cxd > o" > 0. Introducing 
an improved approximation with more poles necessarily brings in non-trivial zeros. 
Can the zeros change dramatically the smoothness of the MHA? So far, all the cal- 
culations made in the literature in the MHA are based on the working assumption 
that the smoothness, beyond the MHA, persists; in other words, one is assuming an 
underlying hypothesis of smoothness of large-A^ QCD which, in full generality, has 
not been proved from first principles. What follows in the next sections is a test of 
this assumption, albeit in a very particular case. 

We can only suggest a possible scenario on how this smoothness assumption may 
be a generic property of large-Ac QCD. It is quite clear that in the case where the 
spectral function is positive, the fact that M{a) is then logarithmically convex |22| 
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provides the required smoothness property; however, in most cases (hke the left -right 
correlation function we are considering here) this property of positivity does not hold. 

In fact, in the case of the left -right correlation function, Witten has proved pT|| , 
under rather general assumptions, that IIlr{Q'^) > for < < oo. This follows 
from the positivity of the measure in the gluonic path integral and Schwartz type 
inequalities of the fermion propagators, Witten's proof, however, is not in general 
applicable to other generalized two-point functions. 

A property of smoothness would follow if one could guarantee that the positions of the 
zeros of the generalized Green's functions are all in the negative half-plane Re z <0 
(i.e. the half-plane which includes the Minkowski axis). This would be the case if 
the polynomial of zeros n^i(-2 + "^j)' ^i^^ ^he polynomial of poles Hdl-^ + '^^^ 
also a stable polynomial ^. 



3. Ulu in the case oian — V — A — V Spectrum 



It is claimed by some of the authors of refs. |12, 13, 14] that the reason why their phe- 
nomenological analysis of the chiral condensates give the same sign for {Oq) and (Os) is 
due to the fact that the hadronic r-decay spectrum is sensitive to the presence of the p', 
while the MHA ignores all higher states beyond the first axial state. Partly motivated by 
this claim ^, we want to analyze here the case, beyond the MHA, where an extra vector 
state V', and therefore one zero a, are also included. Let us collect the relevant equations 
corresponding to this case. 



1. The Correlation Function 

With a spectrum of a pion pole, and V, A, and V' states, the relevant correlation 
function is 

-^nLRiQ') =^ Wlr[z] = A. ^ , (3.1) 

[z + l){z + pA){z + pv) 

where 

Ai^— = ^ = PF and {Oq) = —F^MIM^, = -M^-pfPaPv'- (3.2) 
PAPv My a a 

Because of the positivity of WlrIz] for Rez > 0, the position of the zero has to be 
in the Minkowski axis and, therefore, a > 0. 



^In the mathematical literature this goes under the name of the Routh-Hurwitz theorem. See e.g. ref. [ p3[ . 
We have recently been able to prove that the successive polynomials IljLil-^ + ^^j) for = 1, 2, • ■ • up to 
an arbitrary, but finite A'^, are indeed stable polynomials. The proof, however, lies outside the context of 
this paper and will be published elsewhere. 

^Notice, however, that there are other phenomenological analyses which using the same r-data find 
opposite signs for (Oe) and (Os) [|, 0, 
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2. The Linear Constraint 



Equation ( |2.9D now reduces to the simple relation 



a - {1 + PA + pv) 



1 m 



(3.3) 



M2 (Oe) • 



This is one of the key equations of our discussion, which already provides a semi- 
quantitative argument in favor of the opposite sign option for the condensates (Og) 
and (Og) . The equation states that for (Og) to have the same sign as (Oq), the 
position of the zero has to be far beyond the largest V' -pole: 



Fixing the position of the poles at the values of the observed spectrum (and ignoring 
errors for the purpose of the discussion), one has 



Mv = 0.776 GeV , Ma = 1.230 GeV {pA = 2.5) , My = 1.465 GeV {pv = 3.6) ; 



which means that for the equal sign requirement option to be satisfied one must have 
a > 7.2. In GeV units this corresponds to a mass of 2.1 GeV. Now, as already stated 
at the end of section I, in writing a large- A'^c ansatz for the M^Li?[-2] function, one 
is implicitly assuming an effective cancellation between the extra poles and zeros in 
the complex z-plane which lye beyond a disc of radius sq covering all the poles and 
zeros retained in that approximation. The result a > 7.2 implies that the radius in 
question has to be ^/sq > 2.1 GeV. A priori that seems a good thing because the 
OPE-matching is now applied at > so; i.e. in a more asymptotic region than in 
the case of the MHA ansatz; however, it also implies that there are no further poles 
in the region between Myi ~ 1.5 and the effective mass M^- ~ 2.1 GeV corresponding 
to the zero at a ~ 7.2. This, however, is in contradiction with the observed ai-like 
state at Ma' ^ 1.64 GeV and p-like states at My ^ 1.72 GeV and My" ^ 1.9 GeV 
below Mo- — 2.1 GeV. Alternatively, if one excludes those three states A' , V" and 
V" as all the phenomenological analyses using r-data do in fact, then the position of 
the zero a should be a < -j^ 4.5, implying according to eq. ( |3.4[) , that (Og) and 
(Oq) must have opposite signs, in contradiction with the claims of refs. |l^, |14|. 

3. The Slope Constraint 

This is the relation between Lio and the position of the poles and the zero 



which was already discussed in the previous section. Since Lio is relatively well- 
known phenomenologically it gives a constraint between pp, pA, Pv' ^■ 



a > I + PA + Pv ■ 



(3.4) 



(3.5) 




(3.6) 
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4. The Electromagnetic 7r+ - 7r° Mass Difference 
Recall that 

'i rv 1 TOO 

^l+U = T-^ / dQ^ [-Q^Ulr{Q^)] . (3.7) 
^'^ ^0 Jo 



In the MHA with a spectrum, 



2 I (MHA) 3 a 2 1 1 c 

mt+\ = -—MtfPA logPA- (3i 



In the case of a V-A-V spectrum we find 

2 |(v-A-v)_3a 2 PV ^P^' ~ ^^^P^ ~ ^ + ~ ^-^^(^ ~ ^) 

(3.9) 

which, for a = pv', reduces to the MHA expression. Since the Am = — m^o 
mass difference is dominated by its electromagnetic contribution, we can use its ex- 
perimental value as a further constraint on p^, pv and a. 

5. Reconstruction of the Spectral Function 

In full generality, as shown in the previous section, one can reconstruct the spectral 
function from the knowledge of the zeros and poles and the normalization A]\f. In 
our case, using ( 2.3l|) and ( ^.32 ), this results in 



-ImUiRit) = -lmUv{t) - -ImUAit) , (3.10) 
vr vr TT 



with 



hmUA{t) = F^6it) + FS^- ^ -At-Ml) + -/^hit-so)il + ---), 

vr a {pA-l){pv' - PA) 167r^3 

(3.11) 

and 

ilmnWt) = Fo^^ \ . P^^'l^-'^ ,M t-M^) + -. ^—^ -,Sit-M^,)] 

TT a Upa - l)(/oy' - 1) [pv - ^){PV' - PA) J 

N 2 

+ Y^3^(i-^o)(l + •••), (3.12) 

where the dots in these equations stand for pQCD ag-corrections. Notice that these 
spectral functions, by construction, satisfy the 1st and 2nd Weinberg sum rules. For 
a = Pv they reduce to the spectral functions of the MHA case. 

6. Radiative Widths 

For a meson V in the lowest octet of vector states, the width of the electronic decay 
y — > e^e~ is given by the expression 

Tv^e+e- = ^fvMv ■ (3.13) 
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For an axial-vector A, the width of the decay A n-f, in the chiral hmit, is given 
by PI 



r^^., = ^fl^MA . (3.14) 
The relation between these decay rates and pA, pv and a is as follows: 



a {pA - l)[pv' - pa) 
a [pyi - l){pv' - pa) 



7. The Electromagnetic Pion form Factor 

In Large-A'^c QCD, the electromagnetic form factor of the pion has a particularly 
simple expression 

y ^0 ^ 

where Fy = fyMy and Gy = gyMy are standard couplings of the large-Ac effective 
Lagrangian of narrow states [l^. Requiring that the form factor falls as an inverse 
power of = —t in the deep euclidean, one gets the constraint: 

Y,fvgvM^ = Fi. (3.19) 

V 

On the other hand the slope at the origin of the pion electromagnetic form factor 
determines the Lg coupling jl^ as follows: 

L9 = \Y.fy9y- (3.20) 

V 



The two constraints (3.19) and ( p. 20 ), when restricted to a vr, V, A, V' spectrum, 
become: 



fvgv + fvgv'Pv = pf , (3.21) 

fvgv + fvgv = 2-L9 . (3.22) 

In the chiral limit, the coupling gy is related to the p — > vrvr width as follows 

^P^-- = 9v^^- (3-23) 
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The y'-width, however, (and hence the couphng gy/) is poorly known. We can 
ehminate gy between the two equations (3.21) and ( p. 22 ), which results in a useful 
constraint between Lg, pA, pv'-, and gy (which can be fixed from p — > tttt). 



Matrix Elements of the Qj and Qg Operators 

The four quark operators in question are the ones in eq. (|1.3|). We are interested 



in the evaluation of the matrix elements of these operators between an incoming 
i^T-state and an outgoing 27r-state. 



Evaluation of (Oi) As discussed in ref. @, the {{ttt^) i=2\Q7\K^) matrix element, 
to lowest order in xPT is related to the vev 



as follows 



M, ^ ((vr7r),=2|Q7|i^°> = -^(Oi) 



On the other hand, in ref. |6|, it was also shown that 



{Oi 



(2vr)^ 



(3.24) 
(3.25) 

(3.26) 



where the integral has to be evaluated using the same MS-renormalization pre- 
scription as used for the evaluation of the corresponding Wilson coefficient of 
Qt- One then finds that in large-A^c QCD 0] 



327r2 



M2 



(3.27) 



where = p? exp{l/3 + k) with k depending on the renormalization scheme: 
K = —1/2 in NDR and k = +3/2 in HV. For a ir-V-A-V' spectrum this can be 
written as a function of pA, pv and a in the following way 



3M6 py, 



PA 



a-1 A^ 

{pA-l){py,-l) °^M2 



+PAPV' 



-p\ 



a - py 



1 

log 

{PV - 1)(/5V" - pa) Mlpy> 

A2 - 



a - PA 



log- 



{pA-l){pv' - pa) M^pa 



(3.28) 



For pyi = a this expression reduces to the MHA discussed in ref. 
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Evaluation of (O2) This is a vev which appears in the short-distance behaviour 
of the IIlr{Q'^) function; more precisely 

hm {-Q'ULRiQ')) X Q4 = 47r2^ ^(Oa) + ^{Oi)) + O (^) . (3.29) 

It is related to the evaluation of matrix elements of the Qs operator as follows 

Ms = (M,=2|Q8|i^°) = 4(02) . (3.30) 

As discussed in ref . ||6| , to a good approximation and including next-to-leading 
ag-corrections we have 



and in the case of a ir-V-A-V' spectrum 



X 



^ _ / 25/8 (NDR ) \ a.(/u) 
y 21/8 (HV ) y vr 

(3.31) 



r w Pv \ 2 cr- PA 

{02[n)) = -nT—j-TpF— { pa- 



IQirasip) 0- [ {pA-l){pv' - Pa) 

a — 1 a — Pv 

'PA-, TTT TT - PAPV 



' {PA - l)(py' - 1) {pv - l)(Pv' - pa) 

^ _ ( 25/8 (NDR ) \ a,{p) 
\ 21/8 (HV ) / vr 



(3.32) 



Again, in the limit where pyi = a this expression reduces to the corresponding 
MHA expression discussed in ref. 

9. Duality Constraint 

This can be formulated as the requirement that in the chiral limit, there is no l/Q^- 
term in the OPE of the Adler function defined as 

A{Q^)= rft -lmny(t). (3.33) 

The Adler function is not an order parameter of S^SB and, therefore, it has contri- 
butions from the perturbative continuum. Then, in the case of two explicit V and 
V states plus a continuum spectrum, the requirement in question reads as follows 

c2.f2 , o^2 ,^2 _ 4 / 3iVcas(so) 



2/^M^ + 2/^,M^, = ^^.o^l + ^^^^^^ + ---J . (3.34) 

Using the 1st Weinberg sum rule, this can be written as a simple constraint between 
PF-, PA, Pv, cr and the onset of the continuum sq which, obviously, has to start 
beyond the py-pole; i.e., sq > MpPv'- 

( PV a -PA \ Nc 2 So f S Ncasisp) \ 

''V^-ipA-l)ipv-pA))=^^3KV^8^^^---) ■ ^'-''^ 
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4. Numerical Analyses and Conclusions 



1. The case of a tt - F - ^4 spectrum (MHA) 
(a) Fixing the free parameters pp and pA- 

Confronting the MHA approximation to tlie experimental values of the observables 
introduced in the previous section will allow us to test its consistency and adjust the 
two free parameters pp and pA of this approximation. We use as input the following 
set of experimental data 



= 4.5936 ± 0.0005 MeV , 
Lio = (-5.13 ±0.19) X 10"^ 
Vp^e+e- = (6.77 ± 0.32) X 10"^ MeV , 
r^^^^ = (640 ± 246) X 10"^ MeV , 
Lg = (6.9 ± 0.7) X 10"^ , 
Mp = (775.9 ±0.5) MeV; 



ref. |4| 
ref. H 



ref. 14 1 



ref. 14 1 



ref. [|5| 



ref. HI 



(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



and make the reasonable assumption that these observables follow gaussian proba- 
bility density functions (p.d.f.). In fact, some of these observables, when expressed in 
terms of the MHA parameters depend not only on pp and pA but also on (m^+ ±m^o) 
and/or My. Therefore, it is more appropriate for our purposes to use in our fit pro- 
cedure the dimensionless quantities: 



m„+ + m^o 3 a PA log{pA) 

— Om^ = - — - 



4 vr /9A - 1 



Mp P-"'^'' 

J_r - 



1 f 1 

-PF 1 + — 

4 V PA 



Aira PA 
3 PA - 1 ' 

24pA-l ' 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 



This, however, has the drawback that the three ratios ( [4.7D , ( ^.9] ) and ( 4.10| ) may no 
longer have simple gaussian p.d.f. In order to check this, we have computed their 
p.d.f. ^ and found that for each of them, there is practically no numerical difference 
between the calculated p.d.f. and the gaussian one. This justifies the use of the 
standard statistical regression method to fit our parameters pp and pA with the 
result 

pp = (12.36 ± 0.35) X 10~^ and pA = 1.464 ± 0.004 , (4.12) 



*We do this following the example given in ref. 
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with a Xmin ~ ^-^1 ^ degrees of freedom (dof). The covariance matrix is given by 



The quoted statistical errors of pA and pp have been obtained by evaluating the la 
standard deviation via the solutions of: A.x^{pi) = minj^yj [x^{pi,pj) — Xmin )"■'-• 

We deduce from these results two conclusions: first that the MHA framework is 
statistically relevant and second that the fitted free parameters have small statistical 
errors. Moreover, they obey a multivariate gaussian p.d.f.. For My = (775.9 it 
0.5) MeV, we find Fq = (86.3 ± 1.2) MeV and Ma = (938.7 ± 1.4) MeV, where 
the errors are only the statistical errors of the fit. The corresponding perturbative 
threshold so defined by Eq. (|3.35|) is 1.3 GeV. 

(b) Predictions of the MHA 

Using Eq. ( |2.11| ), we have evaluated the first few condensates of lowest dimension. 
The quoted numbers are given in Table 1 below in the entry MHA (the second line). 
They have been obtained using a Monte-Carlo simulation which takes into account the 
statistical correlation between the two parameters. The second error is an estimate 
of the systematic theoretical error which we have made in the following way: 

i) We give a systematic error to the parameter pp of 0{T p/Mp). Notice that this is the 
parameter which modulates the large-A'^c counting in all our theoretical expressions 
(both for the MHA and the MHA + V). 

ii) We enlarge the experimental error of Mp, which is the quantity modulating the 
dimensions of the calculated observables, by a factor of twenty i.e., 



iii) As already mentioned in the text, the MHA + V framework reduces to the MHA 
one for a = py . This suggests a way of introducing an extra systematic error to the 
MHA results versus the MHA plus an extra pole (R) results, by fixing the a priori 
ignorance one has on the relative position of an extra pole versus an extra zero within 
reasonable limits. We propose to quantify this error as follows 



using the experimental values for pA and pR, {pA = 2.5 and pR = 3.6). Notice that 
this covers the possibility that the extra pole is of the V-type (cr > p/j) or of the 
A-type (cr < PR). 

The three sources of systematic errors are then added in quadrature. 

For the purposes of comparison we also show in the same Table 1 the results of 
the other determinations of the chiral condensates; in particular the values quoted by 




(4.13) 



Mp^ Mp = (776 ± 10) MeV . 



(4.14) 
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Table 1: Numerical Results for the Chiral Condensates 





(Oh) 
xlO^ GeV^ 


(Or) 

xlO^ GeV^ 


(0^n) 
xlO^ GeV^° 


IZ/ 

xlO^ GeV^2 


\>-^i4/ 

xlO^ GeV^^ 


xlO^ GeV^^ 


MHA + V 


-7 90 ± 20 
±1.62 


±11 69 ± 32 
±2.53 


-13.12 ± 0.43 
±3.01 


±13.21 ± 0.62 
±3.24 


-12.54 ± 0.93 
±3.29 


±11 45 ± 1 50 
±3.21 


MHA 


-7.89 ±0.23 
±2.01 


±11.71 ±0.34 
±3.08 


-13.18 ±0.41 
±3.61 


±13.33 ±0.42 
±3.83 


-12.78 ±0.43 
±3.86 


±11.89 ±0.40 
±3.78 


ALEPH 


-7.7 ±0.8 


±11 ± 1 










OPAL 


-6.0 ±0.6 


±7.5 ± 1.3 










Davier et al. 


-6.4 ± 1.6 


±8.7 ± 2.4 










loffe et al.\ll\ 


-6.8 ±2.1 


±7 ±4 










Zyablyuk [15 


-7.1 ± 1.5 


±7.8 ±3.0 


-4.5 ±3.4 








Bijnens et al. 


m 


-3.2 ±2.0 


-12.4 ±9.0 










Cirigliano et 
ALEPH[13| 


2l. 


-4.5 ±0.83 
±0.18 


-5.7 ±3.72 
±0.64 


±48.2 ± 10.2 
±2 


-160 ±26 
±5 


±426 ± 62 
±14 


-1030 ± 140 
±30 


Cirigliano et al. 
0PAL[|l3| 


-5.06 ±0.89 
±0.12 


-3.12 ±3.82 
±0.45 


±38.7 ± 10.6 
±1 


-132 ±27 
±3 


±354 ± 66 
±6 


-850 ± 150 
±20 


Rojo et al.\l^ 


1| 


-4 ±2 




78 ±24 


-260 ± 80 







Cirigliano et al. in ref. |^] ^. Notice that these authors also find an alternance of signs, 
but opposite to our MHA prediction, except for the lowest dimension condensate. 



Using Eqs. ([S^SD and ( |33o| ) we can also evaluate the matrix elements of Qi and Qs 
operators. The results are given in the line MHA of Table 2 below. Again, the first 
error is statistical, the second error is our estimate of the theoretical systematic error 
in the way described above. 

The results in Table 1 and Table 2, corresponding to the MHA, are perfectly consistent 
with those previously obtained in refs. |l^ and |27] using a different treatment of the 
input parameters. 

2. The case ofatt-V-A-V spectrum (MHA±V') 

(a) Fixing the free parameters 

Adding a vector resonance V in the spectrum extends the number of free parameters 
from two to four: pp, pA, pv and a. Furthermore, Eq. ( [4.11 ) becomes now a 



function of gy as explained in Section HI. 7. The way we treat this is by considering 



^See the original reference for a discussion of the two sources of errors. 
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Table 2: Numerical Results for the Q7 and Qg Matrix Elements 





M-j (NDR) 


Mj (HV) 


Mg (NDR) 


Ms (HV) 


MHA + V 


0.12 ±0.00 ±0.01 


0.59 ± 0.01 ± 0.06 


2.00 ± 0.03 ± 0.20 


2.15 ±0.03 ±0.22 


MHA 


0.12 ±0.00 ±0.02 


0.59 ±0.01 ±0.11 


1.99 ± 0.03 ± 0.36 


2.15 ±0.03 ±0.39 


Donoghue et al. \2't 




n 1 ft _|_ n 1 




o 09 ± 67 


2 Afi + n 70 


Bijnens et al. |12| 


0.24 ±0.03 


0.37 ±0.08 


1.2 ±0.8 


1.3 ±0.6 


Cirigliano et al. |^9[ 


0.22 ±0.05 




1.50 ±0.27 




Narison Q 


0.21 ±0.05 




1.4 ±0.35 




RBC H 


0.27 ±0.03 




1.1 ±0.2 




CP-PACS et al. HI 


0.24 ±0.03 




1.0 ±0.2 




Donini et al. p^] 


0.11 ±0.04 


0.18 ±0.06 


0.51 ±0.10 


0.62 ±0.12 


Bhattacharia et al. 


35] 


0.32 ±0.06 




1.2 ±0.2 




SPQCDR H 


0.24 ±0.02 




1.05 ±0.10 





the observable®: 



J_r 



1 cr{pF- 2LqPv,) {pa - l){pv' - 1) 

487rp2, (1 _ p^,)2p^p^p^,(-o- _ 1) 



(4.16) 



as a function of Lg which has an error itself, and is added as an extra parameter in 
our fit. The number of d.o.f. does not change since Lg is also taken as an observable. 
We also impose a criterion of rejection through the ordering: 



PA < PV < (y < Po 



So 



(4.17) 



The first and second inequalities reflect the knowledge that the new state has a higher 
mass than the axial, that it is a l^-like pole and, therefore, its residue contributes 
positively to the W2,ij(z)-function; the third inequality follows from the requirement 
that the perturbative threshold sq defined in Eq. (3.35) already lies beyond the radius 
where the analytic structure of the poles and zeros retained satisfies the leading OPE 
constraint. 

A similar statistical analysis to the one in the previous subsection with a regression 
leads to the following results: 



PF = (12.36 ±0.03) X 10"^ 
PA = 1.466 ±0.003, 
pv = 2.63 ±0.01, 
a = 2.64 ±0.01. 



(4.18) 
(4.19) 
(4.20) 
(4.21) 



^The decay width of p ^ tvtv is Tp 



= (150.4 ± 1.3)MeV 
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The results in Eqs. ( p^ to (|]2l|) correspond to a value: Lg = (6.44 ±0.02) x 10-^ 
with a Xmin ~ 0-^0 "^o^- '^'^^ errors, which are only the statistical errors of the 
fit, were calculated using the same "reduced-x^" procedure as before. We find that 
the parameters pp and pA are statistically stable when compared to those found in 
the MHA case. We also find that py ~ o", which is consistent with the fact that the 
MHA approximation seems to have already the bulk of the full large- A/'c information. 
In other words, adding an extra V'-pole appears to be compensated, at a very good 
approximation, by the position of the nearby zero. 

In order to make numerical predictions for the chiral condensates and the matrix 
elements M7 and Mg we need to know the shape of each "reduced-x^" so as to im- 
plement a Monte-Carlo simulation. Concerning pp and Lg we find parabolic shapes, 
i.e. gaussian behaviours. In the other cases the structure of the "reduced-x'^" are 
slightly more complicated but since, as already seen in the MHA case, the resulting 
statistical errors remain very small as compared to the systematic errors of the theo- 
retical framework, we have finally assumed that all our parameters are gaussian and 
uncorrelated. For My = (775.9 ± 0.5) MeV, we now find Fq = (86.1 ±1.1) MeV, 
Ma = (939.4 ± 1.1) MeV and My = (1258.2 ± 2.5) MeV, where the errors are only 
the statistical errors of the fit, and a perturbative threshold in the ^-channel (or 
A-channel) starting, again, at ^/sq ~ 1.3 GeV. 

(b) Predictions of the MHA+V 

The restriction of (|!8|) to the P = 3 case (MHA±V') reads 

(02n+4) = 2(-l)"+i M« E (1 - -) n 7 ^^/'^ ' (4.22) 

where pi = 1, p2 = Pa and p^ = pv' ■ The resulting values of the chiral condensates 
are given in the first line MHA ± V in Table 1. The matrix elements of Q7 and Qg 
are now obtained using Eqs. (|3l25| ), (|^ and (|3.30D ,(|3^) and the corresponding 



results are given in the first line MHA ± V of Table 2. The systematic errors of the 
MHA ± V predictions have been estimated using the prescriptions i) and ii) already 
described earlier for the MHA predictions. Within errors, the two set of predictions 
from MHA and from MHA ± V are perfectly consistent with each other. 



The shape of the function in Eq. (3.1) with the parameters fixed at the central 



values resulting from the fit in Eqs. ( 4.1 j ) to (|4.2l| ) is shown in Figs. 2a, 2b and 2c. 



Figure 2a shows the shape of the function Wl_r[z] in the euclidean region {z > 0). 
Figures 2b and 2c show the shape of ReWLij[z] (the thick solid red lines) in the 
Minkowski region for —2 < 2; < in Fig. 2b and in the region of the V in Fig. 2c. 
Notice the different scales of the three figures. The delta functions contributing to 
Im Wl/j[2:] are also shown (the thin vertical blue lines V and A in Fig. 2b, and the 
thin vertical blue line V in Fig. 2c.) 

3. Comparison with other Determinations 
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Fig. 2a The predicted shape of the function Wlr[z] in Eq. (3.1) in the euclidean 
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Fig. 2b The predicted shape of the function Wlr[z] in Eq. ( ^.]\ ) in the Minkowski region 
—2 < z < 0. The vertical lines are the delta functions contributing to the imaginary part, 
the continuous solid lines, are the shape of IieWLFt[z] in that region. 



The authors of refs. |12[| and [|13[ give results for quite a few observables. This offers 
the possibiUty of making a comparative study with our predictions. The way we do 
that is by asking the following question: "what are the values of pp, PA, Pv o-''^d a 
which, using the large-Nc parameterization MHA + V given in the text, can fit the 
predictions of these authors?" 

Concerning the work of Bijnens et al. [12|, we have used their predictions (with their 
errors) for (Oq), (Og), M7, Mg and the perturbative threshold (which in the case of 
Aleph they take at sq = 2.53^q }2 GeV) as input data. Using this input, and imposing 



the constraint pA < Pv < cr < po it is possible to find a fit with a x 
and the central values: 



2 

min. 



11.9/(3 dof) 



PF = 0.011,, PA = 1-69, pv'=2.08, a = 4.10, 



(4.23) 
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Fig. 2c The predicted shape of the function Wlr[z] in the V region. 



with 



81.2 MeV and 



Po 



PL 



So = 4.34. 



(4.24) 



The problem with this fit, which reflects the rather bad Xmin. ' is that a < 1 + pA + pv 
in contradiction with the equal sign value for Og and which these authors find. 

It is possible, however, to find a reasonable fit, using as input the predictions of 
Bijnens et al. [^], if one leaves free the position of a. We have found a solution with 
a Xmin ~ 1.2/(3 dof) and parameter values: 



PF = 0-OU7t'oZl (4-25) 

PA = 3Atli (4.26) 

py, = 3.7+};° (4.27) 

a = I5.55tl-ll (4.28) 

Fo = 86.9iy MeV (4.29) 

/Jo = 5.18, (4.30) 



where the errors are only the statistical errors of the fit. The problem of this fit, 
however, is that the position of the zero is far beyond the onset of the pQCD threshold 
{(J ^ Po) and many states, which have not been included in the analysis, can fill this 
gap. 



Concerning the work of Cirigliano et al. [13|, we have used their values for the 
condensates from Oq to Oie, with their errors, as well as their predictions of Mj 
and Mg, with their errors as well, as input values. Like in our previous analy- 
sis, we have first imposed the constraint pA < Pv < Po^^ and a < p^f^'^' for 
So = {1.95, 2.15, •• • ,3.15}GeV^. We have found a solution with, however, a very 
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bad Xmin ~ 134.5/(4 dof). The central values of the resulting parameters are 

/9F = 0.010, p^ = 1.60, pv' = l.78, £7 = 3.48, (4.31) 

with 

Fq = 93.7 MeV and po = 3.59 . (4.32) 

Again, we find that this solution satisfies the relation a < 1+ pA+ Pv in contradiction 
with the equal sign value for Og and Og which the authors find. 

If we relax the constraint on a < pQ it is then possible to find a "better fit" with a 
Xmin ~ 3.26/(4 dof) and central values for the parameters: 

Pir = 0.009, PA = 1-02, py, = 2.03, cr = 5.40, (4.33) 

with 

Fo = 92.39 MeV and po = 3.47 . (4.34) 

The problem with this fit, however, is twofold. On the one hand the position of the 
axial state is too near to the first vector state {pA = 1-02) and also the fact that 
there is still a large gap between the onset of the pQCD continuum (po = 3.47) and 
the position of the zero {a = 5.40). 

We conclude from these analyses that the results of Bijnens et al. [^], as well as 
the results of Cirigliano et al. ||l^, if interpreted within the framework of large- A'c 
QCD, show internal inconsistencies. As recently discussed by Donoghue ||3^, this 
could very well be the reflect of unphysical extrapolations in the hadronic spectral 
function which (implicitly or explicitly) have been made in the phenomenological 



analyses reported in refs. p3, 03, 14 1. 
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